Système de coordonnées ellipsoïdales « confoncales » (notation Jacobienne cohérente avec NIST 
« Hanbook of Mathematical Functions », notation modifiée par rapport à Moon et Spencer) 





Soit le système de coordonnées À,u,v tel que : 
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On peut aussi définir le système en inversant le signe des paramètres A,u,v , puis x et z, comme 
suit : 
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Selon la première donnée et en posant les expressions suivantes, il vient : 
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Les trois iso-surfaces du système de coordonnées ellipsoïdales sont : 
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x y Z 2 X y Z al a. 2 
= a = = -b 
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Lorsque le facteur d'échelle w est égal à 1 il vient : 
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Les trois iso-surfaces du système de coordonnées ellipsoïdales (lorsque y est égal à 1) sont : 
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4 
y =a Ee -e Jen*(a,k)en(B,k)en*(y,k) et 10-30 








2 2 
2 =-0 E-e) a TT 2 
Ce qui donne par exemple comme choix de variable : 


x=@ k sn(&,k)sn(B,k sny, k) 
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Dans ce cas les iso-surfaces s'écrivent: 
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En utilisant les expressions suivantes : 














â =K(k)+iK(k')-a  B=K(k)+ip 

e dn(a,k) = ko 1l 
INE k en(a,k) GE: k en(a, k) 
sal, k)= = L) onl, k) =k' ALL, 





Belo2k(k)] et  7e[0,4K(k)] 
y kz 3% 
sn’(@,k) cn*(X,k) dn (&,k) 
` y 2,2 Le 
sn°(B,k) cn (B,k) an (B,k) 
y 2,2 "a 
sn* (y k) cn’(y,k) dn (7.k) ` 
~); _sn(a,k) 
dn(&,k)= -ik ER) 


ASA p B) 





dn(B,k’) 


On arrive aux expressions du système de coordonnées : 
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Et les iso-surfaces s'expriment finalement comme suit : 
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On peut encore changer légèrement la paramétrisation pour arriver au résultat suivant : 
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On retrouve la forme donnée par Moon et Spencer, par la suite des substitutions suivantes : 
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Avec les iso-surfaces suivantes : 
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Et celle de Morse et Feshbach par les substitutions suivantes : 
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avec les iso-surfaces correspondantes: 














Ce qui donne moyennant une substitution à l'aide des fonctions elliptiques de Jacobi : 





-od =W >") 
dn(a,k) SA -d 1 DS _ Í sn(a,k) y ala, ) 


en(a,k) en(a,k 5 ( bd Gx) 
nla, k dnl B, 
E sadn) 78 dell Sa-i ape)  |2709 >C 04) 


E =bsn(y,k) yb’ -E =ben(y,k) Va -é =adn(y,k) ael0,K(k)] et Bel-Kk(k)K(x)] et y e[0,4K(k)] 











— 








Résultat que l'on a déjà établit précédemment par d'autres substitutions. 


Une dernière forme est également possible à partir de celle de Morse et Feshbach 
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C'est par exemple cette dernière forme que l'on peut utiliser pour définir le système ellipsoïdale de 
dimensions supérieurs. 


Système de coordonnées sphéro-conique (sphero-conal ou sphero-conical en anglais, notation 
cohérente avec Moon et Spencer , NIST) 














2.2 (b: -ï# =) 0<v° <b eu <e et œ facteur d'échelle 
0<r<to 








Cette forme originale donnant les carrés des coordonnées cartésiennes dispense de la nécessité 
d'introduire les signes permettant de définir une bijection sur les six compartiments de l'espace 
cartésien. Par ailleurs on parvient à la forme Jacobienne en posant : 
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On arrive aux notations Jacobienne du NIST Handbook of Mathematical Functions : 


x=0krsn(B,kn(y,k) 0 <r < +œ 
sets B =K(k)+iB 

y=oi rr ren(B,ken(y, k) B < [0,2K()] 

Z=@ Er aoB Alt E) ve [0.4K (K)] 





Les trois iso-surfaces du système de coordonnées sphéro-coniques sont : 
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B =K(k)+iB 
~ 1 = sn(B,k') = cn(B,k') 
,k)= Eh dn\b ,k |= k'——— 
diag Ge POE mga) 
Il vient : 
1 
x =ok „k 
x= o krsn(B,k)sn(y, k) KEEN "us ml ) 
y=0iLrenB,klenr.k) = y= okr MB only à 
V le (B,k') 
z=0 =r dn\B,k dn(y ,k =>. dau dnly E 
k' í Jan ) Z TT Aën nly, ) 
Et les expressions des iso-surfaces sont les suivantes : 
2 2 2 
x DE 32 
r 
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2 Y kz 
= =0 
x tya sn’(B,k') k" cn*(B,k') 
x? y k? z? 
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Autres formes du système sphéro-conique 


Les expressions de paramétrisation du système sphéro-conique que l'on peut rencontrer des formes 
sont très diverses, à tel point que l'on a parfois du mal à vérifier qu'il s'agit bien du même système 
de coordonnées. Prenons par exemple la forme originale et par la suite des substitutions dans les 
fonctions elliptiques de Jacobi, arrivons aux formes suivantes: 


























on(B,k!) | | L z=or sn(K(k') 
z=0r- GE) dn(y,k) dn(K(k')- B,k')= GE 
Be [0,2x(k')] x=or Aë. sn(y,k)  [x=or an(ÿ, DÉI k') 0<r<+to 
3 B= K(k')-Be[-K(k'), K(k) > 4 y=or ol Kenly.) <= y=or enlB,k)cn(7,k') avec B =€ [- K(k), K(k)] 
y € [0,4K(k)] z=or srl, k) dn(y,k) z=or sn(B,k)dn(y, k') y € [0,4K(k')| 


La toute dernière expression (après substitution k et k') est la forme donnée dans l'article de 2010 
R.MENDEZ-FRAGOSO.E.LEY-KOO « Lamé Spheroconal Harmonics in Atoms and Molecules ». 


Pour ces deux dernières formes, les iso-surfaces s'écrivent pour la première forme comme suit : 





























> _ cn(B,k') a enlB,k 2, „2 q 72 
Si snlB,k')= mee R LE er Ce e 
x=or dn(B,k')sn(y,k) m É v) Si , | 
Ss R p oul. sn\p; La snp, X d = A = 
y= or en l)e ,k) > BE te MB) A de Ge) arr) Gr) 
z=0 r sn(B,k')dn(7,k) e k k 2 2 AE) 
dr(B ,k')= => © dn(B,k')= => ege 
dn(B,k') dn\B,k' snè(y,k) cn°(7,k) dn*(y,k) 
Et pour la deuxième forme : 
X+Y +2 =, 
x=or anl Bk) nly.) E F 2 
E Aë lä Al sm (Bk) ` 
z=or sn\p, nly, ) d E vi 2 kez 7 
sn’(y,k') cn’(y,k') dn’ (y,k') 





Partant d'une de ces expressions, introduisons les amplitudes de Jacobi pour arriver à une forme 


dites « trigonométrique », comme suit : 


sn(y,k)= Cos(o) 
Ss cn(y,k)= Sin(p) 
x=0 r dn(B,k')sn(y,k) gp ==-amly,k) dn(y,k)=+/1=k*Cos*(@) 


y=or en(B,k')cn(y,k) 





Et les expressions des iso-surfaces sont les suivantes : 


x+y tz 


x=orJl—-k” Cos*(9)Cos(g) 
=r Sin(8 )Si 
E in(9)Sin(g) 
z =œ r Cos(9) 1-k’Cos’(@) 
O<r<+o = 0<9<x  0<p<2r 





= @ 
Z 
X=or 1-k° Cos*(9)Cos(p) Së y z 
= o r Sin(9)Si 8 j 
y = or Sin(9)Sin(9) + Tk Co (9) Si (9) Cos*(8) 
z =œ r Cos(9) 1=k*Cos*(p) x y kz 








C'est la forme qui est donnée dans l'article de 1973 de 


Cos (eo) Silo) 1-k'Cos*(p) 


l.Lukacs.Ya.A.Smorodinskii e SEPARATION 


OF VARIABLES IN A SPHEROCONICAL COORDINATE SYSTEM AND THE SCHRODINGER EQUATION 
FOR A CASE OF NONCENTRAL FORCES » avec le changement de notation: 3? ¥' & pv 


Par la substitution k et k' depuis la première forme, une modification dans la prise en compte d'une 
amplitude et l’interversion de x et y, on peut également écrire : 


Gë ld x=or Sin(9)Cos(p) 

2 y=or 1-k*Cos*(9)Sin(p) 
z =œ r Cos(9) 1-k" Sin (p) 
0<r<iw 0<9<r  0<p<2r 


x=or dn(B,k)sn(y,k)) 
y=or enlB,k)cn(y,k') @ =am(y,k') 
z=0r sn(B,k)dn(y,k) x & Y 


Et les expressions des iso-surfaces deviennent : 





x =œ r Sin(9)Cos() r? z 
3 {js ) x? k y z 
= or J-FCos (9)si z = 
E "VE 2 Sila). 1-FCos (9) CosT9) 
z =0@ r Cos(9)/1-k" Sin” (o) 2 y pz 








Cos'(p) Silo) 1-k" Sin'(9) 


Ou bien encore les substitutions suivantes (en prenant des angles complémentaires dans les 
amplitudes de Jacobi) : 





sn(y,k)= Sin(o) 
S en(y,k)= Cos(p) = ` GO Qu fac 
x= or n E = x=aryl-k” Sin (9)Sin(p) 
3 @ = amly,k) dn(y,k)= J1- Sin (o) y =œ r Cos(9)Cos(p) 
y=orcnlp,k)en(y,k = \> = < l => 
Ge 9 =am(B,k') ` | sn(B,k)=Sin(9) 2=0 r SONT ES (p) 
Z=orsn Jk!) dnly,k) ~ 
enlB,k')= Cos(s) O<rsto 0<9<2n 0<p<2r 
dn(B,k')= 1K” Sin"(9) 


Et les expressions des iso-surfaces sont les suivantes : 


REY Es 
2 
x=or 1=k" Sin’(9)Sin(g) Se S z 
de 
d- 7 sin*(9) ` Cos*(9) _ Sin*(9) 
2 2 


x y Ee 


Sin*(p) R Cos*() 4 (1=kSin*(9) 





y =œ r Cos(9)Cos(p) > 
z=or Sin(9) 1-k*Sin" (9) 


=0 











Une dernière forme que j'ai trouvé dans la littérature est due à R.D.Spence, mais ce n'est 
probablement pas la dernière possible, s'obtient à l'aide de la forme Jacobienne, comme suit : 


cn? (8, k) 


> 
>» 





a 
ana) 877 dn°(y,k 





sn (ÿ,a) SC sn*(y,.k') 








k +Ç" =dn(B,k) JE +n? =dn(y, k) 





DN x=00, dëi Clg all 





et 
B? 


>a’ +n 





Les iso-surfaces de cette paramétrisation sont les suivantes : 














r 
rm, dÉ +E LÉI =n? Peya 
B 242 
r r 
J> x? Eë z? BR 
? Y 7 2 2 2 DT 
r EC, $ E 
mg, do -6" ja +n’ SZ A a 
0<r <+œ ao’ =0/=07=1 B'-n' n o +n 


C'est la forme que l'on trouve dans l'article de 1959 R.D.Spence « Angular Momentum in Sphero- 


Conal Coordinates » 


Système de coordonnées paraboloïdales générales 





La forme originale est donnée par les expressions suivantes : 


























SCC fee 

a i À = -0<V<b</À<a<u<œ 

J ag, [Oe „ç ub -bX =b) 0<b<a et o facteur d'échelle 
o ` a-b d) b-a ko y 

zi o, =0, =l 

<=-(u+2+v-a-b) 

o 2 





Notez que dans les définitions usuelles, outre que x et y sont exprimés au carré (ce qui introduit un 
signe sur x et y), la variable v n'est pas négative, or pour qui veux regarder de plus près, sur une 
iso-surface u=Cste, si v est strictement positif et comme À reste toujours borné, alors 
nécessairement z est borné ce qui n'est pas possible car l'iso-surface est une nappe de paraboloïde 
infini, donc avec z non-borné. La correction de cette erreur systématiquement reproduite est 
simple, il suffit de considérer v comme pouvant être infiniment négatif. Et dans ce cas toutes les 
zones de l'espace peuvent être couvertes avec chacune des familles d'iso-surfaces orthogonales. 


Les iso-surfaces s'écrivent comme suit : 





a-u b-u 
2 2 
Kory -0[2-1] 
a-A b-A o 
2 2 
x e y SERA 
a-v b-v o 


En posant (voir P.Morse.H.Feshbach-Methods Of Theoretical Physics Part |, page 664) : 























= 2 12 
pe ee Stee ae Aie ae), od 
a Va cn’(y,k) cn’(a,k) 
2 12 
=ü = p2 Sn" (v. k) ÀA-a=-ak" sn7(B.k = _ __ ak" ` 
DEER Steg a=-ak"sn*(B,k') v-a PE 
k" *(a,k 








= q sn°(B,k'}sn°(7,k) 
cn*(a,k)en*(y,k) 
k 


D sn(B,k')sn(y,k) 
o d 
y? 5 sn* (a, k)cn"(B, ') 

@ 


cn(a,k}en(y,k) 








2 
o 
= SC _ g sn(a,k)en(B,k') 
en*(a,k)en*(y, k) o en(a,k}en(y,k) 
z d| sn(a,k) dn'(B,k) sn(y,k) z dl sn(a,k) dn'(B,k) sn°(y,k) 
=—|- + + =—|- + + 
oa 2\ en’(ak) k? = eny,k)) |œ 2 en(a,k) = k* ois 











Comme on peut rassembler en un seul terme le facteur d'échelle w et d, il vient : 


R+k?=1 kel 
a e |-K(k), K(t)] 








pelox(r] EH 
y e K(k), K(k) Ü _ sn'(a,k) _ dn*(B,K') sn’(y,k) 
2\ cn’(a,k) ke n° (y,k) 


De même en posant : 














* = (a-b) Sinh(a)Sin*(B Cosh’ (y) 
Bet = Sinh°(a ) et (a) e 
= ae Z = (a-b) Cosh (e)Cos*(B)sink?(y) 
= Cos’(B) TL = Sin’ (8) > o 
ab E z (a Hl sua) — Sinh’ (y )+ Cos’ (b)+ 2 
<= Cosh’(y) ëch De a-b 
< 8 = 2. lH cosiaa)-cost(27)+Coslop))+ Ë 


Il vient alors en posant le facteur d'échelle w au profit d'un nouveau facteur d'échelle c comme suit 
et en translatant d'une constante fixe l'axe z : 





2 
@' (a-b) =4c x? = 4c’ Sinh*(œ)Sin*(B)Cosh*(y) > x = 2cSinh(œ)Sin( B)Cosh(y) 
pe => ï y" =4c’Cosh’(a)Cos’(B)Sinh’(y)=> y = 2cCosh(a )Cos(B)Sinh 
5 y B > y B)Sinhly 
a _ 
z- cath z =5 (Cosh(2a)- Cosh(27)+Cos(2p)) 
2a-b 


Enfin par une permutation de x et z, il vient : 


x= 2cSinh(a )sin(B \Cosh(y) = 5 (Cosh(2a a)- Cosh(2y)+ Cos (28 ) 


y= 2cCosh(a)Cos(B )Sinh(y) > y = 2cCosh(a)Cos(B}Sinh(y) 
Z = "(Cosh(2a)- Cosh(2y)+ Cos(2B)) SE 
2 0<a<œ O0Sy<w -r<y<r 


C'est la forme que l'on retrouve aussi bien dans l'article de 1967 de F.M.Arscott e XXI—The 
Whittaker-Hill Equation and the Wave Equation in Paraboloidal Co-ordinates » que dans la thèse 
de 1968 de K.M.Urwin « PARABOLOIDAL WAVE FUNCTIONS ». 


Système de coordonnées sphéro-coniques à N+1-dimensions 


Soit un système de N+1 valeurs réelles ao a, ..., Ane, et soit les coordonnées sphéro-coniques 
1S1,...,Sn telles que : 


dd a, ER; et  r,5,.,5,€R; 


X >0,x,>0„..xy>0 x, 





>Vx; ER x;=60;r 





Iso-variété dimension N 











A trois dimensions on retrouve bien les coordonnées sphéro-coniques déja définies moyennant les 
changements de notations suivantes : 








ae. h-e 1) ; , (b? ul =v) K et ` œ facteur d'échelle 


Bec) ` > Pe-e) 0 <r <+œ0 











z =@r (u -c*\v* =c?) _ 2 le-w fe =v?) 
ele -b) cc? -b°) 
x? =p (s, als, -&) 
(a, "Ou Xa, = a) 
XX, Mäh 2x a =0 a=b a =c > 2 (s,-a,Xs,-a,) 


5 A Dix =r 
V >s. Uu >s, =l ( 








Bonne nouvelle, l'équation de Laplace à N dimensions est simplement séparable dans le système 
sphéro-conique : 








i=N 1 i=N 1 i=N-1 EE 
(s =a, w."(s,)+ _ (s ) + Y As, lw (s )=0 
J. LAIN E JO NET j Vi i? j jN j 
i=0 i=0 $j — 4; i=0 
PEN A 
we 
op | D As, 
" ' i=0 = 
ow; (s,)+ 2 Teg (s,)+ i=N w(s,)=0 
i=0 
i (s, d) 





A trois dimensions ce sont des équations de Heun et même de Lamé pour ainsi dire. Les équations 
différentielles sont des équations fuchsiennes à N+2 singularités do, a, ..., Ane et l'infini. 


Réalisons maintenant la projection stéréo-graphique de la sphère S"* de rayon 1 sur l'espace R" 
selon son pôle Nord : 


Pour i>0 x =Â x; =- 


1-x, T a,-a,) 


Les iso-variétés de dimensions N-1, sont alors définies à l'aide des iso-variétés du système sphéro- 
conique à N dimensions en utilisant la projection stéréographique, comme suit : 

















LN 
Go 
i=N 3 3 i=N j Sr J É x 
ales [D x; = |X, +> x avec x; =y x =-=>+@& x, =(1=x,)x'; 
i=0 il 7 ) 1-x, 
d a; 
j=0,i+j 
i=N 
1 
e eh oc iN x; ed 2 iN lex 
= ! — ' ee = | 12 _ 0 

> X, =1-(1=x,) Yx, SIE EEN D X = a eS ee et EE 

F p 2 T s — 

Ge S Ss Lag z 0 

i=l i=] 
p | ) 
x'<-1 
i=N 2 2 i=N 2 12 t i=N 12 

x x l-x,} x = 4 x 

j e I... N) i = _ `= dre pe = + i_=0 








Nouvelle bonne nouvelle, le système de coordonnées est R-séparable cette fois-ci. Ici je vais 
reprendre la construction réalisée par H.Volkmer et H.S.Cohl dans leur article de 2013 : 
« Separation of variables in an asymmetric cyclidic coordinate system ». Prenons tout d'abord la 
projection stéréo-graphique à N+1-dimensions, et considérons les coordonnées r et x'; comme une 
paramétrisation de l'espace cartésien à N+1 dimensions, comme suit : 


(XX pees xy JES" = dE = bar d 


P:S” \(1,0,.,0)>R" tq Batz, ae Gee) 





=N 
| 2, | 
PR > S" \(1,0,...,0) tq P(x px Je —— 


(oe SKI o 
: Zo 





Les coordonnées curvilignes (r,x') présentent comme facteur d'échelle (de Lamé) les valeurs 
suivantes : 


h=1 h=- = A(x, „....x'„ )=2rh 


be -1,2x', = j i F? Ps 


i=N > E 
2 

X x +1 

i=l 








(xX Xy ) = r 
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=0 
E ax. ry, a, | h? ôx, 








Valeur du Laplacien qui est la suivante: x y= 1 |ô [am ZLS ô DA" òV 
(ark) | ar ôr ] F òx', ôx' 


Si l'on suppose comme le font H.Volkmer et H.S.Cohl que la fonction U est homogène de degré a 
alors la fonction V peut être décomposée comme suit : 


Ob x tx „tx ls U(x x ,.. xy) => V(r, Xsan PN W(x" res x) 


a-2 | i=N 
S Apah =" 52 le AME = < "wen IN W(x' Ee oX 2 
i=l i 


4h" x", 





En introduisant la fonction u comme suit, alors l'expression du Laplacien se simplifie encore : 














OW _ p Ou Oh Ou 
> =_ = pL S =h — -2x ne 
is k ox! Co +up Or, ' E X; up 
Ox", 
3 Ô hr? ow Ô kär: 2 ôu -26x', uhh L = 
Ôx', Ôx', ` e. Go : 
ie =E -2x\ hP (2B+N- T — 2B uhh 44px? (BEN -Iuh 
X ; 
Si 2B+N-2=0 
N ô [owr0W] 31 ôu : AH 
SE => h"? =h? -2Buh? +2NB x'; uh? 
f 2 Z 2 ax’, p a 
B+N-2=—-1 
i=N ôW(x' N i=N x. (=N. Ny 
== e: SS sx d h? +o -2NBuh? +2NB) x' uh? 
FI x" i Ô i=l 





1 
i=l Or, ox’; i=l i=l 


i=N i=N ' ' N i=N ai Ni i=N 92 
diets VL dp ont) =n? S'S" op unt re S2 NW 2)uh? 
h ' Or Ka 
De cette dernière expression, on tire la valeur du Laplacien : 


re? Ny i=N EN u 3 Na 
Ay V =| h? EEN +4a(a+N-I)h? u 


i=| OX 








a-2 i=N ail re? i=N a2 
= by Var] 22 4 (4a(a+N-1)+N(N-2)Wu =+ S YZ 4(N+20)N- E 
4h? i=l Ox", Ape i=l Ox", 


Autrement dit pour deux valeurs données du paramètre a, si la fonction U est harmonique dans le 
plan cartésien R"* alors cette fonction s'exprime à l'aide d'une fonction également harmonique 
dans l'espace projeté R" et inversement : 


i=N =l 
U(x,, Xi Xy ) = rh SO ike) KH 
i=l 








EN 92 
Unes) DE =06 a=-7 ou a=1-7 
i=0 OX; 
i=N a2 
A u =Y 2 “=0 





Dans ce cas on voit donc clairement que le Laplacien est R-Séparable dans l'espace projeté par 
projection stéréo-graphique de la sphère de rayon unité, avec un système quelconque de 
paramétrisation orthogonales sur cette sphère (ici je reprend le système sphéro-conique, mais cela 
aurais pu en être un autre) : 





On remarque également que dans le cas N=2, soit la projection stéréo-graphique de la sphère 3D 
sur le plan, le Laplacien reste séparable car la fonction de séparation est égale à 1. 


Le système sphéro-conique à N+1 dimensions, sur l'iso-variété de la sphère de rayon unité, une fois 
projeté par projection stéréo-graphique devient R-Séparable. La fonction U peut s'écrire dans les 
coordonnées sphéro-conique à N+1 dimensions comme une fonction homogène en r de degré a. 
Ce dernier doit avoir une valeur fixe pour respecter les conditions de R-Séparabilité lors de la 
projection stéréo-graphique : 


u(x x J=w(s)w,(s,) ef Si a=-— ou Ge 


l= 


X'y X JE T” ms}.w,(s,) = A„V =0 





EI 
Fonction de séparation R(x XX AIS x a i) 


Ce système de coordonnées est appelé N+2 cyclides car les équations différentielles séparées 
comporte N+2 singularités. Du fait de l'homogénéité de la fonction radiale, cela entraîne une 
contrainte sur la première constante de séparation : 


nb watch ab y (r)=06 ole Aa Na +44, =0 
r P 








ER 2(N-1)- N° 
, Seu [TT 7777 16 
= 
4 ach PEN , a = NW-2) 
2 2 16 


Si N=2Sa=-1>4=-< ou _a=0—1,=0 


Si N= EVE XT AA: ou nn Ae 
2 16 2 16 


Dans l'espace à trois dimensions c'est l'expression de cyclides asymétriques. Si l'on applique cette 
construction à l'espace à 4 dimensions, en coordonnées sphéro-coniques, puis que l'on réalise une 
projection stéréo-graphique sur R°, on arrive à une paramétrisation de l'espace R°. Soit : 


2 2 2 2 
a <s, Sa <s, Sa, Ss, Sa, 0 =0, =0, =0, =l 
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2 
+4 ea = +4 =0 
S; - Ay Zä Sja, ZO 


En prenant comme convention ao=0,a;=1,a2=a et a3=b, on arrive à la paramétrisation N°12 donnée 
par W. Miller « Symmetry and separation of variables » page 210. La solution séparée de 
l'équation de Laplace s'écrit dans ces coordonnées comme suit : 


AU = E CH _, ux, p, z) = V (Sra (sahus (s) 


2 
ôx ` Ô a? (lie ay te" 
S=S,,85,8; W=W,W,,W, A,=-—ou A = ee constantes de séparation 


Gen SL kuke aen ut 
Leen) ES uk Ze erst 





Prenons do=0,a1=1,a2=2 et a3=3, voici trois exemples d'iso-surfaces, d'abord s1=(a5+10a:)/11 : 


| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
| 
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Puis S2=(a1+5a2)/6: 





Puis 53=(a2+Q3)/2: 





La version à deux dimensions est la suivante : 














x, =0, E al 
\ (a =h a, -&) i x, 
fx = 
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2,42 _1Ÿ 2 2 
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Iso — courbes cartésiennes asymétriques 


La solution séparée de l'équation de Laplace s'écrit dans ces coordonnées comme suit : 


OU OU 
AU at ag =0 U(x, y)= m (s )w,(s,) 
ox” ôy 
1 
S=8,9  W=W,W AE e ou Ass As Ad, constantes de séparation 
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On reconnaît une équation de Lamé dans les deux équations différentielles. 


Système de coordonnées ellipsoidales à N-dimensions 


Soit un système de N valeurs réelles ordonnées a, ..., Ay, et soit les coordonnées ellipsoidales 
S1...,S\ telles que : 


+ + 2 
4 „..,ay ER; et S,.,sy ER, ef o; =1 














2 SH j=l = 
0<a, EA EA x = >VxeR x=0;r 
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i=1 $j 4; 

A trois dimensions on retrouve bien les coordonnées ellipsoïdales déjà définies. Encore une bonne 
nouvelle, l'équation de Laplace à N dimensions est simplement séparable dans le système 
ellipsoïdales, avec un système d'équations différentielles fuchsiennes du second degré identiques à 


celles du système des coordonnées sphéro-coniques à N+1 dimensions (l'indexation des constantes 
de séparation est simplement inversée) : 


An constantes de séparation s;=>s ww 








i=l i=l Sj =d; 





